Thermodynamics of a Bose-Einstein Condensate with Weak Disorder 
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We consider the thermodynamics of a homogeneous superfluid dilute Bose gas in the presence of 
weak quenched disorder. FoUowing the zero-temperature approach of Huang and Meng, we diago- 
nahze the Hamiltonian of a dilute Bose gas in an external random delta-correlated potential by means 
of a Bogoliubov transformation. We extend this approach to finite temperature by combining the 
Popov and the many-body T— matrix approximations. This approach permits us to include the 
quasi-particle interactions within this temperature range. We derive the disorder-induced shifts 
of the Bose-Einstein critical temperature and of the temperature for the onset of superfiuidity by 
approaching the transition points from below, i.e., from the superfluid phase. Our results lead to a 
phase diagram consistent with that of the finite-temperature theory of Lopatin and Vinokur which 
was based on the replica method, and in which the transition points were approached from above. 



I. INTRODUCTION 

An interacting ultracold dilute Bose gas in a weak ran- 
dom external potential, which is homogeneous in the 
mean, represents an interesting model for studying the 
relation between Bose-Einstein condensation and super- 
fluidity and has been the subject of various theoretical 
investigations in the last few years P, H, d, |3) B [l] ■ Two 
different methods have been used for performing the av- 
erage over the impurity scatterers. In Refs. P, [2, 0] the 
average of the grand potential over the disorder is taken 
perturbatively in the strength of the disorder after di- 
agonalizing the Hamiltonian by means of a Bogoliubov 
transformation. Alternatively, in Refs. the avera- 

ging is implemented by the replica method. In Ref. 
the replica symmetric solution of the model is found by 
a systematic diagrammatic Beliaev-Popov perturbation 
theory for the dilute superfluid gas in the presence of dis- 
order. At zero temperature the two different approaches 
give equivalent results. At flnite temperatures, however, 
the approach based on the Bogoliubov transformation 
becomes unsatisfactory since it necessarily neglects im- 
portant quasi-particles correlations P, The replica 
trick, on the other hand, has limitations of it's own by 
making the mathematical and also the physical descrip- 
tion less transparent. For these reasons it would clearly 
be desirable to develop an alternative theory avoiding 
these shortcomings. This is the goal of this paper, where 
we show how the perturbative approach of Refs. [U, 0, Q 
can be extended to include the leading effects of the scat- 
tering between quasi-particles at finite temperature. As 
a result, some discrepancies between the two different 
methods are resolved. 

The paper is organized as follows. In order to make the 
paper self-contained, we briefly rederive in Section HIl the 
thermodynamic potential and the equation of state for 
the case of a vanishing spatial correlation length of the 
disorder potential In the limit of zero temperature, 
we also determine the high-order Beliaev corrections to 
the chemical potential [4| . In Section IIIII we give the 



derivation of the superfluid component of the system in 
the presence of a random potential. The disorder-induced 
corrections of the velocity of sound are there obtained 
from hydrodynamic equations for the superfluid compo- 
nent of the system. In Section [IV] we extend the theory 
to finite temperatures within the mean-field Popov ap- 
proximation. In Section |V] the thermodynamics and the 
phase diagram are investigated by means of the many- 
body T— matrix approximation. In particular, we use the 
theory to calculate the shift of the Bose-Einstein critical 
temperature and of the temperature for the loss of the 
superfluidity when approaching the critical points from 
below. Comments and conclusions remain for Sect ion IVll 



II. BOGOLIUBOV'S THEORY 

We consider the effects of an external random fleld on 
the thermodynamics of a dilute Bose gas. The random 
field is assumed to have a probability distribution P[U] 
normalized to one when averaged over all disorder con- 
figurations, that is J d[U]P[U] = 1. The average over 
the disorder fields is defined as 



d[U]>P[U], 



(1) 



and for the disorder potential we assume the ensemble 
averages 



(C/(x)) = 0, 
(C/(x)[/(x')) i?(x- 



x') 



(2) 



We assume that any time scale of the disorder potential is 
frozen, i.e. very long in comparison with the thermody- 
namic time scale. This so-called quenched-disorder limit 
has the consequence that the disorder average must be 
taken after the thermodynamic average over the grand- 
canonical ensemble. Therefore, the total average of an 
observable reads 



d[u]p[u]{om,.. 



(3) 
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where (0(^'))gr indicates the grand-canonical average. 
At very low temperatures, when the wavelength of the 
atoms in the gas is much larger than the range of the 
impurity scatterers responsible for the random potential 
Q, one may consider the limit of a (5— correlated type of 
disorder i? (x — x') = _Ro 5{x — x'). The constant Rq is 
then related to the concentration of the impurities and 
to the s— wave scattering length of the random scatterers 

0- .. . 

At equilibrium, the grand-canonical partition function 
of a Bose gas in a disordered medium is given as a func- 
tional of the random external potential U (x) 



j_S[r,t,U] } , (4) 



where the functional integral is performed over c-number 
fields il!*{yi,T) and ipi'^:''') periodic in imaginary time 
over h(3 = h/ksT. The Euclidean action S is given by 



temperatures, the de Broglie wavelength Ath of the atoms 
is much larger than the range of the interaction poten- 
tial such that ro/Ath ^ 1. Therefore, only s-wave scat- 
tering is relevant in the gas and we can neglect the mo- 
mentum dependence of the interaction potential, setting 
Vq ^ Vq. Breaking the gauge-symmetry of the action by 
introducing the decomposition -(/"(x, t) = y/rio + tp' (x, t) 
and expanding the resulting expression up to quadratic 
order in ip' and U, the effective action becomes 

S^^'> [a, a*] = ~hf3fmaV + ^hPnlVoV + hf]noUo (8) 



k.i 



+ 



E 

k.n 



-k.-r 



^k.n^ — k,— r; 



(5) 



OT 



2m 



Ai + C/(x) 



/ dr dx '!/'*(x, t) 

+ iydx' V*(x',T)y(x-x')^(x',r)|v(x,r) , 

where /i is the chemical potential and ^(x — x') the 
atomic repulsive interaction potential. The disorder av- 
erage of the thermodynamic potential f2 = — (InZgr) /P 
is obtained from 



(InZg,.) 



(6) 



The average of Eq. ^ is highly non-trivial because the 
disorder average is nonlinear in Zg^ due to the logarithm. 
In this paper we calculate this average by following the 
method of Huang and Meng in Ref. [l| which is based on 
a canonical Bogoliubov transformation Q . 

Expanding the fields in Fourier modes in 
a finite volume V according to ip^x^r) — 

(l/n/3y)^/^^k.n"k,ne*(''-''^""^) with the Matsub- 
ara frequencies w„ = 2TTn/hP, and the corresponding 
complex conjugate expression for '0*(x, t), we write the 
action in momentum space as 

S[a*,a] =^{-ihujn + eu- fJ.)al „au,n (7) 

k,n 

+ y ^ <^k,n^k-k'ak',n(5n,0 



k.k'.n 



2h(3V 



q '*k+q,»i+m'^k' -q,»i' -rn "^k' ,n' Ok.n 



k,k',q 

n^n' ,711 



In this equation we have introduced the free particle dis- 
persion — ft'^k^/2m and the Fourier transform of the 
interaction potential Vq = / dx y(x)e~*'''^. At very low 



where the prime denotes that k = is excluded from 
the sum. The condensate density no remains to be de- 
termined by minimizing the thermodynamic potential. 
Note that we have also neglected as of higher order the 

term l/^I]k,k',n"k,ra t^k-k' ak',n^n,o, which requires the 
assumption of weak disorder [T]. The effective action is 
diagonalized by a Bogoliubov transformation 



flk.ri = MkOk,! 

at „ = utal 



Wka-k,-n - ^k 

z* 



^k,n - "k"k,ri - vt<^~K~n - Z^k7 (9) 

where the coherence factors Uk, Wk and the complex num- 
ber Zk can be taken real positive by appropriately choos- 
ing the phase of the complex fields. In such a case we 
have 

Ck - /i + 2?ioVb' 



2 1 
^ 2 



1 



2 1 

^ 2 



-1 



Ck 



fj. + 2noVo 



t/k 



z±k = [no/V)^ — — (wk - Wk) 
n.S2k 



(10) 



where |ukp — l^^kP = 1, Wk^k > 0, and where the Bogoliu- 
bov spectrum is given by 



ftfik = Y (ek - M + 2noVo) - (rioVo) 
After the diagonalization the action reads 



(11) 



1. 



a, a 



+ 



E 

k.n 



+ -^2^ [f^k - (ek - /i 

k 

/3^^'|jy-^|2Mek - A 



2710^0) 



V 



(12) 
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By performing the functional integral of Eq. (|4]) and av- 
eraging over the disorder, the thermodynamic potential 
of Eq. dHl) becomes 



(O) = - ^,naV + \nlVoV + ^In (l - 
+ [/if^k - (eic - M + 2noVb)] 

k 

En (<^k - + noVo) 



(13) 



The first two terms on the r.h.s. represent the mean-field 
result while the remaining three terms describe thermal 
and quantum fluctuations. From the thermodynamic re- 
lation n = — (l/V) d{n) /dfi we have for the particle- 
density in the grand-canonical ensemble 



n — no + n' + tir, 
with the depletion due to the normal interaction 



(14) 



-y' 



k 

ek 



" ek - M + 2noVo 



HQ,], 



and the disorder-induced depletion 



-e' 



noR( 



(fk - M + noVp) 



(15) 



(16) 



k 

Here 7V(ftf]k) = (e'^'^"'' - 1)"! is the Bose distribution 
function. The chemical potential can be eliminated from 
this expression by minimizing , given in Eq. , with 
respect to the condensate density uq [cf. Eq. pTjl below] . 
In the mean-field approximation this yields fi ~ noVo, 
which corresponds to the Hugenholtz-Pines relation for 
that approximation. At this minimum, the thermody- 
namical potential of Eq. becomes (with ng ~ f^/Vo) 



(0) 



(17) 



[^f^k - Sk - noVo] - ^ 



RqUo 



ek + 2noVo ' 



with hfli^ ~ Ek (ck + 2noVb)- Analogously, in the equa- 
tion of state of Eq. ^4]) , we get 



no: 



noa 



V 



ek + TtpVo 



N{hnk), (18) 



and 



-e' 



(ek + 2710^0)'' 



(19) 



Performing analytically the integration in Eq. ((T9)) . the 
number of condensate particles depleted by the disorder 
becomes [H 



riR = Ro 



no 

8 7r3/2 V a 



(20) 



where we have made the replacement Vq — > , that 
is, we have eliminated the unknown bare potential Vo 
for the two-body scattering matrix T^^ = Anh^a/m pro- 
portional to the s— wave scattering length a. According 
to the limitations of the Bogoliubov approximation, the 
theory is valid under the conditions nR,n' <^ no — n. 
The constraint n' <^ n implies the diluteness condition 
n^/^a <^ 1. The condition tir <C n is equivalent to the 
inequality R'o = m'^Ro/Sn'^^^h'^ (na)^^^ <C 1. Associating 
to the strength Ro of the disorder perturbation a length 
scale defined as c? = (27rfi2/m) /Ro, the condition on 



R'o can be rewritten as 2 (n^/^ay ~ dn^^''^ /^/t: ^ 1. In 
the dilute limit n^^^a <C 1 this requires dn^/'^ ^ 1, i.e. 
the length scale associated with the interaction energy 
due to the impurity potential must be much larger than 
the interparticle distance n"^/^. The condition i?Q <C 1 
can also be reexpressed and elucidated by introducing 
the healing length of the condensate wave function. The 
latter is defined as ^heai = 1/V87ma. According to the 
theory of superfluidity in BEC, the inverse of the heal- 
ing length characterizes the upper boundary of the mo- 
menta in the phononic spectrum of the fluid. At this 
wavelength, the energy of the excitations is of the or- 
der of ftr^k — M- The condition i?g <C 1 is equivalent 
to V^n^hcai ^ d. Therefore, the theory is valid when 
the energy of the excitations, induced by the impurity 
scattering, is far below the value ft/^hcai ^^^^ marks the 
crossover from the collective phononic excitations to the 
single-particle excitations. 

The lowest-order Hugenholtz-Pines condition fj, = 
noVn, neglects the effects of quasi-particle interactions 
[lOl . as well as the scattering between the quasi- 
particles and the impurities [1]. Nevertheless, the beyond 
mean-field Beliaev corrections to the leading order result 
fi = noVo depend, both in the normal and in the dis- 
order interactions, only on two-body collisions and can 
be calculated in the framework of the Bogoliubov the- 
ory ^l3|. Minimizing the thermodynamic potential of 
Eq. with respect to the condensate density, we have 



^noVo + Vo-Y, 



2ek — 2fi + inoVo 



2ek -2fi + 3noVo - 2hnk 



2?ifik 
ek - /i -I- 2noVo 



(21) 



'2noVo (ck - + "-o^o) (2ek - 2^ + SnoVo) 
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In order to get the corrections to the mean-field result we 
substitute in the right-hand-side of equation Eq. (PT|) the 
zero-loop result fi — noVo to obtain the next order correc- 
tion to the relation between the chemical potential and 
the condensate density. At T = 0, the Bose distribution 
N{hn-k) can be neglected, and we have 



(22) 



k 



£k - 2noVb 
' (ek + 2noVaY 



Subtracting the ultraviolet-divergent contribution 

^^Jol4?(l/^)X]k l/2ek in the third term and eliminat- 
ing in the above expression the bare potential Vq for the 
two-body scattering matrix defined by the Lippmann- 
Schwinger equation T^^ = Vo~ Vq{IIV)Y,^^ ^T^^ , we 
find 



^l =noT2^ I 1 



40 / riQa? 



3 V TT 



T^^nR. (23) 



Using the equation for the number of particles in 
Eqs. (Ull), mi) and (HH), we have that Eq. ^ can be 
rewritten as 



^ = nT2^ I 1 + r: W 1 + ^ (24) 



32 na 



3 



3 V TT 



with 



<5m 



(ek + 2,^oT2S) 



(25) 



in agreement with the zero-temperature result of Ref . [3] ■ 
Note that the beyond mean-field correction given by i5/iR, 
still contains an ultraviolet divergency, which is not re- 
lated to the interaction and was thus not yet removed by 
the renormalization of the latter. Rather its origin lies in 
the fact that we have considered a (5— correlated random 
potential. In second order perturbation theory, it changes 
the original chemical potential /i = /iba 
dom potential, to the new value ^ = /iba 



without the ran- 



Thus, after averaging, we must make the change fi — 

+ Ro{i/V)J2k which removes the diver- 

gency. After renormalizing in this way, we have 



(26) 



which corresponds to a shift for the macroscopic com- 
pressibility 9((5/iR,) /dn — Sx/xo = 3rtR,/n with respect 
to the mean-field value xo = T^^n/m of the Bogoliubov 
theory for a clean system. Using the thermodynamic re- 
lation fj, = dF/dN, we can calculate the free energy F 



from Eq. 1^. At T 
E and we find [Ij, i^] 



this coincides with the energy 



{3_ 

V 



1 



128 



na 

TT 



8i?n 



(27) 



III. SUPERFLUID COMPONENT 

In a Bose-Einstein condensate, random impurities con- 
stitute a source of incoherent scattering which tends to 
localize the condensate. As shown by Huang and Meng 
fH the formation of local condensates in the minima of 
the random potential reduces the superfluid component 
of the fluid even at zero temperature, where, in the ab- 
sence of disorder, the whole fluid would be superfluid [31 ■ 
For the sake of completeness and in order to propose a 
simpler derivation, we rederive in this section the super- 
fluid component in the presence of weak disorder [l], 0, 0| 
by extending the Bogoliubov diagonalization method to 
a moving system. The superfluid density rig is related to 
the total density by the relation 

ns=n - Un, (28) 

where n„ is the density of the normal component of the 
fluid. In order to calculate n„ , we consider the action of 
Eq. ([7]) when the gas is in motion. The moving reference 
system is related to the laboratory system by a Galilean 
transformation x' = x + ui and t' = t. The flelds in 
the moving system have a relative velocity with re- 
spect to those in the inertial frame. They are related by 
the transformation ^*'{x',t') = e-™^^^/'^** (x, i) and 
*'(x',i') = e™^=^/''«'(x,i). Therefore, in the new ref- 
erence system, the action in Eq. ([7]) becomes 



S[a*,a] = ^ [-ihwn + /ik (u - v^) -f ek - Meff] a^n^k,™ 

k,n 

— ^ ak,n^k~k'ak^n^n,0 (29) 
I 

y ^ ^ '^k+q,n+m'^k'-q,n'-m'^k',n'Q^k,n7 



^ k,k' 
1 1 



2hf3V 



k,k',q 

n,n'va 



where the new chemical potential is defined as /Xofi = 
fi + mwVs — TOv^ /2. In the broken symmetry regime, the 
action in the new reference frame, up to the quadratic 
order in the fluctuations flelds, reads 

5(2) [a, a*] = -n/3/icfrnoU -I- ^hf3nlVoV + hPnoUo 



-E 

k,n 



[-ihwn + fik (u - Vg) -I- ek - /Ueft + 2^0 Uo] ak_„ak,„ 

jL / 

nQ\2 \ ^ 



VJ ^ 

k,n 



I I 



(30) 
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The latter action can again be diagonalized by the Bo- 
gohubov transformation described in Eqs. (P|)- (fTT|) . This 
is achieved by replacing — s- ficS in the definition of Uk, 
Vk and SrJk, and by defining the new shift variable zj-k 
as 



z±k = {no/Vy 



± fik (u — Vs 



(uk - Wk)^ • (31) 



Performing the diagonalization, the functional integra- 
tion and the average over the disorder, we obtain the 
averaged thermodynamic potential 



+ [ftf^k - (ck - Moff + 2noVo)] 



^ noRo 



(ck - Mcff + noVo) 



(32) 



Expanding for small ftk (u — v^) to second order, we have 



(n) ~ - ^,,snoV + InlVoV + ^ ^In (l - 
^ k ^ 

+ [fif^k - (ek - Moff + 2noVb)] 



k 



(ck - Mcff + ?io"^o) 



(33) 



2 [e-zsiOk _ 1] 



RqUq 



2ek 



(u- 



where the term linear in u — Vs vanishes as a conse- 
quence of the symmetry hQ]^ — fi,57_k of the Bogoli- 
ubov spectrum of Eq. (|lip . Moreover, minimizing this 
thermodynamic potential with respect to no, we obtain 
t he zero-loop re sult /icff = noVo which gives HQ.^ ~ 
^/ekJek'+'2JI^ . At this minimum, the momentum of 
the system can be calculated from the thermodynamic 
relation 



d{n (T,F,/Xeff (u),u)) 



We find 



p =mVnVs H — y 



+Ro m h k 



2,.2 4ek 

Kit 



(u- v^), 



(34) 



(35) 



where we have used the thermodynamical relation n = 
— {l/V)d{Q) /dfx and the identity ducs/du = mv^ 



Therefore, we can conclude that the density of the nor- 
mal part of the fluid moving with velocity u is given by 



1 2 e'^''"^ 



(36) 



Note that in two dimensions the derivation of the nor- 
mal component due to the disorder is analogous, but the 
factor 1/3 in the formula for the thermodynamic poten- 
tial in Eq. (j36|l is replaced by a factor 1/2. After the 
integration in two-dimensions, we get 



Rqht? 



(37) 



in agreement with Refs. [1, 

The depletion of the superfluid density due to the dis- 
order affects also the propagation of an external distur- 
bance through the system, because the collective motion 
of the superfluid component is "hampered" by the com- 
ponent of the condensate localized in the minima of the 
disorder potential. In order to see this effect, we now 
calculate the corrections to the velocity of sound induced 
by the disorder at zero temperature. Let us assume that 
for weak disorder the dynamics of the superfluid compo- 
nent of the gas can be described by the phenomenological 
two-fluid hydrodynamic equations 



d 

—n + V {vgUs + v„n„) = 



= 0, 



(38) 



where /i is the chemical potential given in Eq. (j24p . How- 
ever, in that expression, we shall neglect the second-order 
Beliaev term due to the normal interactions in compar- 
ison to the corrections due to the disorder, and we shall 
focus on the latter. Then Eqs. ([38|) represent a Landau 
"two-fluid" model for the superfluid condensate and the 
localized non-uniform condensate. What is somewhat 
unusual is that, in this case, the normal component in- 
duced by the disorder is at zero temperature and does 
not carry any entropy of the system. Furthermore, we 
assume that for low frequency excitations, only the su- 
perfluid component can react to the probe, while the lo- 
calized normal component remains stationary. Such a 
situation is familiar from the propagation of the fourth 
sound in ''He [l^ which is expressed by the condition 
v„ = 0. If we restrict ourselves to the linear regime, 
we can write n{t) = n + Sn{t) and /i = fJ-o + Sfi with 
(5/i = [dfi/ dn) 5n. Then, Eqs. (|38p give the equation of 
motion 



m—6n- 



V 



on 



= 0. 



(39) 



From Eqs. ^ and we have that (dfi/dn) = 

rp2B _|_ ^ji-^/jiy Moreover, from the result for the su- 
perfluid density in Eq. we have Ug = n {1 — An^/Sn). 
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Therefore, Eq. ([39|) can be put into the form 



d 



^5n - c^V'^Sn ^ 0, 



(40) 



which exhibits a phonon dispersion Hw = cq. Within 
this direct approach, the sound velocity is found to be 
~ Cq (1 + 5nR/3n) in agreement with Refs. 0, Hi- 
Note that Cq — T^^n/m is the mean- field value of the 
Bogoliubov theory for the clean system. The derivation 
of the sound mode we gave here is rather general and suf- 
ficiently simple to be generalized in order to calculate the 
effects of the disorder on the frequencies of the collective 
modes in trapped gases [l^l- 



Popov theory is equivalent to replacing in the contribu- 
tion riQVb/2 to the pressure ^l/V the bare interaction Vq 
by the renormalized T^^-matrix and to adding the con- 
tributions 2nonT^^ + n^T^^ After these steps the 
averaged thermodynamic potential can be rewritten as 



(n) = - fj.noV + -nlT^^V + 2nohT'^^ + h^T^^ 



1 



IV. POPOV'S THEORY 

At finite temperature the interactions of the thermal 
component of the gas are described by the contributions 
to the action beyond the quadratic order given in Eq. ([5]). 
In this section we extend the Bogoliubov approach of 
Section [IT] including these fluctuations according to the 
scheme of the Popov theory [18] which is designed for the 
temperature domain fc^T > /i. In that approximation 
the cubic and quartic contributions are taken as 



and 



5(3) [a, a* 



no 



h(3V 



(41) 



(42) 



k,n 



where the temperature dependent total depletion n — 
gr) is still to be determined. For our present ther- 
modynamic considerations, we neglect the cubic terms 
following [l^ and include the quartic term. The cubic 
term in Eq. (|4ip only contributes in second and higher 
orders of Vq and is taken into account by the introduction 
of the T^^— matrix below. The new action is still diago- 
nalizcd by the same Bogoliubov transformation Eqs. ([9])- 
PT|) but with the difference that the chemical potential 
/i has now to be replaced everywhere by the new variable 



fi' = fi- 2hVo, 



(43) 



and that the condensate density uq becomes strongly 
temperature dependent. Performing again the functional 
integral and the disorder average, the thermodynamic po- 
tential of Eq. pS]) reads 



(n) - -^lnoV + InlVoV + Y 4ln (l - 

k 

+-Y [ft^^k - (ck - m' + 2noVb)] 

k 

(ek - + noVo) 



(44) 



(fk - /i^ + npT^^) 



(45) 



The equilibrium condition at fixed temperature is found 
by minimizing the thermodynamic potential of Eq. (j45p 
with respect to uq. Using the modified Hugenholtz-Pines 
relation /i' = tiqT^^, which fixes n according to Eq. 
as n = [(/x/T^-^) - nQ]/2, we obtain instead of Eq. ([21]) 



=nQT +7ioT I 

2ek + noT^^ 



Ck 



1 '^R, 



y^"" (ek + 27ioT2S)2 



(46) 



where is given by the same expression as in Eq. (jl9p . 
In the "high" temperature limit fcsT 3> uqT^^ , the main 
contribution to the momentum integral containing the 
Bose distribution comes from the region et < uqT^^ and 
the Bose distribution can be approximated as N{x) ~ 
ksT/x. Therefore, Eq. (|46p can be rewritten as 



/ rp2B / 1 40 npa^ 



1 



ek 



^?^°(ek + 2noT2S)^ 



+ T^^UR. (47) 



In the limit of zero disorder Rq = this result reduces, of 
course, to that of Popov [T8|. The new equation of state 
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1 d{n) _ 

V dn ^ 



no 



8 
3 



no a 



Rotio 



(48) 



This latter equation has the same form as the equation of 
state of the Bogohubov theory as given in Eq. (|14p with 
(fTS]) and (fT9| , but the domain of vaUdity and the details 
of the temperature dependence are, of course, different. 
In the "high" temperature region ksT > noT^^ where 
the Popov theory applies, we can neglect the quantum de- 
pletion of the zcro-tcmperature theory, and the thermal 

depletion n' = (1/T^) Ek'^(^^k) [(ck + noT^^) /hQk] 
can be simplified as 



Therefore, Eq. can be rewritten as 



no 



TO ) ^^+-H- 



(49) 



(50) 



The curve for the critical temperature Tc in the Popov 
theory can be obtained by putting no (T) equal to zero in 
Eq. ([50)1 . We note that the contribution in Eq. ([501) due to 
the disorder vanishes when approaching the critical point 
because, according to Eq. ([^ . we have nR, ex ^JnQ{T). 
Therefore, we can argue that even the presence of a ran- 
dom potential of the kind under consideration here, the 
Popov approximation does not shift the value of the crit- 
ical temperature away from that of an ideal Bose gas T^. 
In the absence of disorder, the Popov approximation of a 
dilute Bose gas describes a first-order phase transition: at 
the critical temperature the condensate density ng (T) 
exhibits a discontinuous jump to a finite value. The lat- 
ter can be calculated analytically [20| as a function of the 
scattering length a. It can be shown, that the effect of 
the disorder is to suppress this discontinuity and that the 
jump vanishes for R'^ larger than some value R'q. 

Nevertheless, the transition to superfluidity occurs at 
the temperature Tg where the superfluid density vanishes, 
which means, where the equation 



V'^ 3 

k 



g"-R, (51) 



is satisfied. At "high" temperatures such that fc^Ts ^ 
noT^^, the latter equation can be approximated as 



TP 



2/ rj.2B\h m^keT, 



(52) 



Therefore, Eqs. ()50p and ([52)1 represent two coupled equa- 
tion for no and Tg. In the lowest order in the strength 
of disorder and interaction we can approximate Eq. ([50]) 



>(3/2) 



. Sub- 



1 - (T/T^) 
and neglecting there the second 



with the ideal gas result no 
stituting this in Eq. ([S? 

term on the right-hand side, in the limit R'q ^ [na ^ 
we obtain an equation for the critical temperature T^ as 
a function of the total density and the strengths of the 
interaction and the disorder 



8 7r3/2 ?l4 



-R'o- (53) 



Solving for Tg we find for the critical temperature where 
superfluidity disappears when coming from lower tem- 
peratures 



T,/T,"~l- (32/27) i?;,' 



(54) 



This result coincides with Eq. (32) of Ref. [1]. Here it has 
been derived assuming that near the critical temperature 
Ts the interaction between the bosonic particles can be 
described by the temperature independent T^^— matrix. 
This is consistent only if the critical temperature Ts is 
not too close to the transition temperature T°. The rea- 
son is that near the Bosc-Einstein condensation point the 
quasi-particle interactions acquire a strong temperature 
dependence 0, |2l| . As pointed out in Ref. 4^ , the 

consistency condition is satisfied when R'q ^ (na'^)^^^. 

In the opposite regime, when i?Q ^ (na"^)^^^, we will 
show in the next section that we have Ts where Tc 

is the Bose-Einstein condensation temperature including 
disorder effects which depends linearly on Rq. In that 
case, temperature effects on the particle-particle scatter- 
ing cannot be neglected and the two-body T^^— matrix 
of the Popov theory has to be replaced by the many-body 
T— matrix, which we shall do in the next section. 

The failure of the Popov approximation in the regime 

Rq ^ [na^Y^^ ti^s a specific physical reason. The Popov 
approximation, as considered here up to now, neglects 
the Hartree-Fock contribution by which the presence of 
the impurities affects the scattering among two thermal 
particles. This can be seen, for example, by calculating 
the relation between the chemical potential /i and the to- 
tal density n and comparing it with the result of Lopatin 
and Vinokur in Ref. [4]. Using Eqs. ((37])-((301), we find 



^ = T^^n - T^^- in')-' 



(55) 



mkgT 



1 



Ek 



(ek + 2fi') 



We observe that in comparison with the theory of Ref. Q 
we seem to miss in Eq. ([55]) the Hartree-Fock term 
T^^ Roui^kBT / Air'^hP . A corresponding contribution 
seems also to be absent in the equation of state in Eq. (|5C 
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and in the normal density component of Eq. ((52)) . How- 
ever, since we are in the regime R'q 3> (na"^) the latter 
contribution is in fact negligible and one obtains the value 
of Tg as indicated in Eq. ([54)1 . In order to access also the 

regime the next section, we will extend 

the Popov theory to the many-body T— matrix approx- 
imation. In this way we will find results similar to the 
one of Lopatin and Vinokur but within a gapless theory 
approaching the critical point from below. 



V. MANY-BODY T-MATRIX 

In the vicinity of the critical temperature of Bose- 
Einstein condensation, the interactions between the 
quasi-particles are strongly renormalized by temperature 
effects and vanish at the transition point for this rea- 
son [2l|. For a clean system, Bijlsma and Stoof Q 
have shown that the quasi-particles interaction is well 
described by the many-body T-matrix T*^^ . In the pres- 
ence of weak disorder, the many-body T-matrix continues 
to obey, in lowest order, the same formal Bethe-Salpeter 
equation as in the absence of disorder [7] 



{1 + N+ + N^) 



2 9 



z — Ml-L- — Ml^ z + + Ml- 



9 9 



9 9 



Z -\- hSl-y — hsi- Z — hfi^ + Ml 

X T*^^(k",k',K;z 



(56) 



where N+ = N{hVt+) and = N{hVt^). The plus 
sign denotes the momentum argument K/2 4- k", and 
similarly the minus sign denotes the argument K/2 — k". 
In the low-temperature domain, where BEC-experiments 
with cold atoms are always realized, the momentum and 
energy dependence of the T*^-^-matrix can be neglected. 
This even applies to the domain ksTc ^ ksT ^ fi of suf- 
ficiently high temperatures where the Popov approxima- 
tion applies, to which we sometimes refer as the "high"- 
temperature domain in this context. Then we have 



T*^^(0, 0, 0; 0) = Vo- VbT*'''(0, 0, 0; 0) 



MB/ 



(57) 



dk 

(2^ 



2nr2L 



4 {hn^) 



[1 + 2Nihn^)] 



The integral on the right hand side contains an infrared 
divergency [l^l caused by the well-known fact that the 
finite-temperature theory does not properly account for 
the dynamics of the phase fluctuations in the infrared 
limit p^ . However, in the "high" temperature limit 

na Ajjj ~ 1 the infrared divergent term must, in fact, be 
dropped. This follows because at "high" temperatures 
the Bogoliubov spectrum ftfJk deviates from — fJ.' only 
for a very small interval of momenta around zero, where 



Ek ^ n-oT"^^ < nT^^. Therefore this phononic part of 
the spectrum merely makes an asymptotically small con- 
tribution to all thermodynamic quantities. Therefore, we 
can set Wk = 1 and Wk = in Eq. With this, and 

after having eliminated the bare potential Vq by means of 
the Lippmann-Schwinger equation for the T^^— matrix, 
we obtain 



T^^'^(0, 0, 0; 0) = T2^(0, 0; -2//) - T^^{0, 0; -2/i') 



(27r)3 hQk 
Solving for T^^'^(0, 0, 0; 0) we find 0, H 

T2S(0,0;-2^') 



T^^^(0,0,0;0) = 



1 + T2B(o, 0; -2fi') J 



(59) 



The many-body T-matrix approximation in a dilute ul- 
tracold Bose gas can now be obtained from the formulas 
of the Popov theory as developed in the previous sec- 
tion just by replacing (20| there the two-body T-matrix 
T2^(0,0;-2/x') = T^^ by the temperature dependent 
many-body T-matrix T*^^(0, 0, 0; 0). The excitation 
spectrum is thereby changed to 



[ek - /-i + ftSii(k, w„)] - [hSi2(k, ujn)] , 



(60) 



where the self-energies are given by 

;iSi2 = TioT^'^ (0,0,0; 0) , 
hY.li = 2nT*^'^(O,O,O;0). 



(61) 



The condition ^' = ^-ftSn = noT^^^(0, 0, 0; 0) ensures 
the spectrum to be gapless. 

The Bose factor in the many-body T— matrix of 
Eq. (|59|) leads to a temperature-dependent scattering 
length, defined as 



AT = mT^^^ (0, 0, 0; 0) /AttH^ 



(62) 



For temperatures not too close to the Bose-Einstein 
critical temperature, where the usual mean-field Popov 
theory of the previous section is valid, the denomina- 
tor in Eq. (|59p represents a negligible correction. In 
that case, the many-body T— matrix indeed reduces 
to the two-body temperature-independent T^^— matrix 
rp2_B _ Anh^a/m. However, approaching the critical 
region near the Bose-Einstein transition, the scattering 
between quasi-particles now becomes strongly tempera- 
ture dependent 0, ■ A more detailed insight can be 
gained by considering the "high" temperature expansion 
of Eq. ([55)1 poj . In that case, the latter reduces to a 
quadratic equation which can be solved analytically. Us- 
ing Eq. (j62p we find for the temperature dependent scat- 
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Figure 1: (a) Temperature dependent scattering length as a 
function of the rescaled temperature T /T^ calculated numer- 
ically from Eq. (|59|) (solid line) and from the analytic expres- 
sion of Eq. (|63p (dashed line). The dotted line shows the 
asymptotic limit described by Eq. (|64|l valid in the Popov re- 
gion for temperatures not too close to T^. (b) Temperature 
dependent scattering length as a function of the rescaled tem- 
perature r/r° near T° . The lower line shows the numerical 
result from Eq (|59|l compared to the asymptotic expression 
given in Eq. (|65p described by the upper line. In both pictures 
the gas parameter of the Bogoliubov theory has been chosen 
such as (na^)^/^ ~ 0.01. 



tering length 



T 



(63) 



where g 



2V^/C(3/2)3 (na^) 



,1/6 



In Figure [T] the 

curve resulting from the "high" temperature expansion 
of Eq. (j63[) is compared with the curve obtained from 
Eqs. ([59|) and ([62| . For temperatures not too close to 
the critical point, we have (qT/T^) < 4(no/n) and the 
curve of Eq. (|63p can be approximated by 



ax — a 



C(3/2)^ 



,1/6 T 



(64) 



while for temperatures just below the critical temper- 
ature the scattering length between quasi-particles 
becomes a universal function of the density and temper- 
ature given by 



ax 



C(3/2)3 „g 



An 



,4/3 ' 



(65) 



which vanishes at the critical temperature T° 0, [2l[ with 
the same power as the condensate density. The two dif- 
ferent asymptotic regimes of the expression in Eq 
are also shown in Figure [1] Note that Figure 1 (a) 



eludes also the zero-temperature limit where the Popov 
approximation as well as the many-body T— matrix are 
infrared divergent. However, this is an asymptotic small 
region because the condition of validity of the Popov ap- 
proximation defined by naX^^ ~ 1 is equivalent to impose 
the condition T/T^ - [C(3/2)^ /ATT]{na^Y'^ . 




Figure 2: The lower and the upper dashed lines describe the 
condensate density no and the scattering length ot obtained 
by solving self-consistently Eqs. (|63p and (|66p in absence of 
disorder, i.e., when R = Q. In this case, both quantities 
vanish aX T = and no shift of the critical temperature is 
induced The gas parameter of the Bogoliubov theory 

has been chosen the same as in Figure [1] The solid lines show 
the solution for finite disorder. We see that the temperature 
at which they vanish is shifted with respect to the critical 
temperature of the ideal gas. In order to make clearly visible 
the effect we have considered in the figure the case R'q = 0.5. 
However, for such a value the condensate depletion due to the 
disorder at zero temperature is already ~ 0.4% of the total 
density. Note also that, in contrast with Popov theory, for 
each value of B!q the curve for no obtained from the many- 
body T— matrix approximation exhibits a second-order phase 
transition. 



By applying these results to the equation of state in 
Eq. (|5D|) the latter can be rewritten as 



n = riQ+n 



T 
TO 



no- 



Anfi^aT 



TO2 ksT 




8 7r3/2 h'^ ■ 



(66) 



For a clean system neither the two-body nor the many- 
body T-matrix theory induce a shift of the critical tem- 
perature [7]. In contrast to this, we find that in the pres- 
ence of delta correlated disorder the effects described by 
the many-body T-matrix induce such a shift of Tc- This 
is illustrated in Figure [21 where the solutions of the cou- 
pled equations ((63|) and ((66)) for ax and the condensate 
density no are shown as function of the temperature in 
absence and in presence of disorder. The origin of the 
shift becomes evident when considering the limit T 
in Eq. (j66p . If we insert in that equation for the temper- 
ature dependence of the condensate density that of an 
ideal Bose gas, the term due to the disorder does not 
vanish at as it would in the simple Popov theory. 
This follows because the temperature dependent scatter- 
ing length flT goes to zero with the same exponent as the 
condensate density no, and both effects cancel. Thus we 
get 



n = n (r/r°) ' + Rom^kBTj^/STT^ff 



(67) 



10 



which is solved by 



^ (1 - ry) , 



(68) 



where 77 = Rom^ksT^ /Uir^ffn = 

27r/[3C(3/2)^/^dni/3] < 1. This value differs from 
the Hartree-Fock result found by Lopatin and Vinokur 
in Ref. [1] by a factor of 1/2. Note, however, that in 
contrast to Ref. Q, our gapless T*^ ■^-matrix theory 
approaches the critical point from below. In fact, the 
theoretical descriptions on the two sides of the critical 
points are remarkably different as a consequence of 
the difference between the two phases on both sides. 
The result of Lopatin and Vinokur has been recently 
confirmed by Zobay in Ref. @ by means of a one-loop 
Wilson renormalization group calculation. This latter 
method approaches the critical point from above as well 
but takes into account critical fluctuations which are 
non-perturbative in the interaction. The numerical solu- 
tion of the renormalization group equations shows that 
the critical fluctuations lead only to small corrections 
with respect to the result of Ref. ; 4] . 

The many-body T— matrix approximation is also use- 
ful to understand the dependence on the disorder of the 
superfluid transition Tg. Equation ([52|) can be rewritten 



TP 



Eq. 

at T 




(69) 



together with Eqs. and evaluated 

Ts constitute a set of closed equations for Tg, 
no\T=Ts ^'^^ 0'T\T=Ts- "^^^ self-consistent solution of the 
three coupled equations is shown in Figure O When 

i?o > (na^)^^^ we have that Tg lies well below in 
the domain, where the many-body T— matrix is given 
essentially by the temperature independent T^^-matrix. 
In that limit is given by Eq. (jM)) of the Popov theory. 

However, when i?o < (na'^)"'^^^, one expects Ts ~ and 
the temperature effects on the quasi-particles scattering 
become important. In that case, we can use the expan- 
sion of the many-body T— matrix near the critical point 
given in Eq. From Eq. to the lowest order 

in the disorder and in the normal interaction this leads 
to Ts ~ T° (1 — 477/3), which is smaller than the critical 
temperature Tc given in Eq. ([55)1 due to the factor 4/3. 

Therefore, we find that the presence of disorder shifts 
the critical temperature for the superfluid transition be- 
low the line of the Bose-Einstein condensation tempera- 
ture. The shift increases monotonically as a function of 

the disorder strength i?o. For R'^ < (na^Y^^ the shift is 
linear in Rq. At larger values of the disorder strength, 
3^l/6 



where Rq 3> 



no") it becomes quadratic in Rq. The 
crossover between the two regimes is described accurately 
in Figure[31 Note that, in that figure, the results obtained 




Figure 3: The solid line indicates the critical temperature 
Ts of the normal to superfluid transition as a function of the 
dimensionless disorder parameter R'q obtained by solving self- 
consistently Eq. ((63l), Eq. IMJ, and Eq. ((Ml) for (na^)^/^ = 
0.01. The upper dashed line is obtained by approximating 
Eq. (|63|l with the asymptotic solution of Eq. (|65|) which is 
valid when Ts — Tc for R'g < (na^)^'^^. In contrast, the 
lower dashed line describes the regime R'q > (na^)^^^ where 



Ts <^ Tc and Eq. (|63|) can be approximated with Eq. (|64|I . 
The comparison between these three curves describes clearly 
how the temperature effects in the many-body T*^^ —matrix 
determine the crossover from the quadratic regime to the lin- 
ear regime for the dependence of Ts/Tc on the disorder cou- 
pling constant Rq. In addition, the short-dashed line shows 
the effects of neglecting the self-consistency in the equation of 
state H66p . where the condensate density is approximated with 

the ideal gas result no = n[l~(Ts/T°f^% For R'o > (na^)^^^ 
the short-dashed line coincides with the analytical result of 
the Popov theory given by Eq. (|54|l . 



in the region Rq ^ 1 have to be considered as an extrapo- 
lation to the more complicate regime of strong disorder 

MM- 

Moreover, the superfluid density can be calculated 
from the relation Us — n ~ Un with the normal den- 
sity given by the right-hand side of Eq. evaluated at 
temperatures T < Ts- In contrast with the Bogoliubov 
approach of Ref. [H , our theory includes important finite- 
temperature correlations between quasi-particles and we 
find that the superfluid density decreases monotonically 
as function of the temperature. This result is in agree- 
ment with the diagrammatic theory based on the replica 
method developed in Ref. Q. 

Here, we have not discussed the damping of the sound 
due to the impurity scattering. This damping has so far 
been studied only at zero temperature in Refs. 0, 
The many-body T— matrix theory developed here, ne- 
glects the finite lifetime of the quasi-particles and cannot 
describe damping phenomena for this reason. 
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VI. CONCLUSIONS AND OUTLOOK 

In conclusion, we have extended in this paper the per- 
turbative approach developed by Huang and Meng in 
Ref. [H for a homogeneous superfluid dilute Bose gas 
in the presence of weak disorder. We have shown that 
such an extension to finite temperatures is achieved via 
a suitable combination of the mean-field Popov approx- 
imation and the many-body T— matrix approximations. 
In particular this allows to make contact with results of 
second-order perturbation theory at finite temperature 
developed in Ref. [3] by means of the replica method. 
E.g. the shifts of the two different critical temperatures 
for the appearance of a condensate density and of a su- 
perfluid density have here been calculated when coming 
from the low-temperature side. 

The theory could have other interesting applications. 



Using the notion of the quasi-condensate [Tl| , Andersen 
et al. [T9I have shown that the many-body T— matrix 
theory can describe the thermodynamics of a clean two- 
dimensional Bose gas near the Kosterlitz-Thouless super- 
fluid transition, and it seems plausible that this continues 
to be true also for weak disorder. Therefore, the approach 
presented here could be useful in order to study the 
long-standing problem concerning the disorder-induced 
shift on the superfluid transition temperature in a two- 
dimensional Bose gas. 
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